Reconciling cyanobacterial fixed-nitrogen 
distributions and transport experiments with 
quantitative modelling 

Aidan I Brown and Andrew D Rutenberg 

Department of Physics and Atmospheric Science, Dalhousie University, Hahfax, Nova 
Scotia, Canada B3H 4R2 

E-mail: aidan@dal.ca, andrew.rutenberg@dal.ca 

Abstract. Filamentous cyanobacteria growing in media with insufficient fixed 
nitrogen differentiate some cells into heterocysts, which fix nitrogen for the remaining 
vegetative cells. Transport studies have shown both periplasmic and cytoplasmic 
connections between cells that could transport fixed-nitrogen along the filament. Two 
experiments have imaged fixed-nitrogen distributions along filaments. In 1974, Wolk et 
al found a peaked concentration of fixed- nitrogen at heterocysts using autoradiographic 
techniques. In contrast, in 2007, Popa et al used nanoSIMS to show large dips at the 
location of heterocysts, with a variable but approximately level distribution between 
them. With an integrated model of fixed-nitrogen transport and cell growth, we recover 
the results of both Wolk et al and of Popa et al using the same model parameters. 
To do this, we account for immobile incorporated fixed-nitrogen and for the differing 
durations of labeled nitrogen fixation that occurred in the two experiments. The 
variations seen by Popa et al are consistent with the effects of cell-by-cell variations 
of growth rates, and mask diffusive gradients. We are unable to rule out a significant 
amount of periplasmic fN transport. 
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1. Introduction 

Filamentous cyanobacteria are model multicellular organisms that form unbranched 
filaments of clonal cells. In conditions of low exogenous fixed-nitrogen, terminally- 
differentiated heterocyst cells develop along the filament separated by clusters of 
photosynthetic vegetative cells [HE]. Heterocysts fix dinitrogen [3] and provide the fixed- 
nitrogen (fN) to vegetative cells where it accommodates ongoing growth. Vegetative 
cells cannot themselves fix nitrogen due to oxygen produced by photosynthesis [1]. 
Because fN transport is required for ongoing growth, and because local fN deprivation 
is a key early signal for heterocyst development, it is important to understand how 
fN is transported within the cyanohacterial filament. Since fN is a metabolite many 
convenient genetic visualization techniques are not directly applicable. Nevertheless, 
radio-isotopes of nitrogen were used to show that fN is rapidly incorporated into 
glutamine (molecular mass 146 Da) within heterocysts [5l[6]. Fixed-nitrogen is likely 
transported to vegetative cells as glutamine or other amino acids [1, 4]. Various imaging 
approaches have indicated cytoplasm-to-cytoplasm connections between adjacent cells 
and a shared contiguous periplasmic space, either of which might transport fN between 
cells [71. 

Early indications for direct cytoplasmic connections between adjacent vegetative 
cells within cyanohacterial filaments were electron-microscopy images of small ( ~ 50A 
diameter) holes in the septal cell membrane called microplasmodesmata [3 El [9]- 
Larger pores have also been observed connecting heterocysts to adjoining vegetative 
cells [To]. Recently, direct cytoplasmic transport of the small exogenous fiuorophore 
calcein (molecular mass 623 Da) was observed in Anabaena sp. strain PCC7120 
(hereafter PCC7120) using fluorescence recovery after photobleaching (FRAP) 
After photobleaching of a cell, fluorescence recovered on a timescale of tens of seconds 

— allowing transport constants to be estimated. Since calcein is not native to 
cyanobacteria, it is thought that direct cytoplasmic transport is nonspecific and could 
be used to transport fN. The protein SepJ (FraG) appears to play a key role in 
direct cytoplasmic transport. SepJ includes a coiled-coil domain that prevents filament 
fragmentation and a permease domain required for diazotrophic growth p2] • SepJ also 
appears to be needed for microplasmodesmata formation [13J. There are no indications 
of ATPase activity of SepJ, or any other energy requirements for transport between cells 

— indicating that direct cytoplasmic transport of small molecules such as fN between 
cells is likely to be by passive diffusion. 

In addition to direct cell-to-cell connections, electron microscopy images indicate 
that the periplasmic space of (Gram negative) cyanohacterial filaments is contiguous 
(see Flores et al [H]). Supporting this, transport of periplasmic GFP (molecular 
mass ~ 27 kDa) has been recently observed in FRAP studies [I5], though this is 
somewhat controversial [2],[71[l6]. Nevertheless, the outer membrane has been shown to 
provide a permeability barrier to fN compounds [17] . Furthermore, knocking out amino- 
acid permeases (transporters) that act between the periplasm and cytoplasm leads to 



Reconciling cyanohacterial metabolite and transport experiments with modelling 3 

impaired diazotrophic growth flEl UHl ED]- Together these imply that there may be 
some amount of periplasmic iN present during diazotrophic growth, which would then 
passively diffuse between adjacent cells via the contiguous periplasmic space. However, 
it is not clear how large diffusive fluxes are in the periplasm compared to those between 
cytoplasmic compartments. 

There have been two studies of nitrogen transport along cyanobacterial filaments 
[3 [21], but neither had the transverse spatial resolution to distinguish cytoplasm 
from periplasm. The first study was by Wolk et al in 1974 [3] (hereafter simply 
"Wolk"). Diazotrophically growing filaments of Anabaena cylindrica were provided 
with dinitrogen gas composed of ^^N (half-life of less than 10 minutes) for 2 minutes 
before being imaged. Radioactive decay tracks from ^^N were used to approximate the 
distribution of labeled fN in the filament. Wolk found that the fN was peaked at the 
location of heterocysts — consistent with diffusive transport away from a heterocystous 
source. 

A more recent study by Popa et al in 2007 [2T| (hereafter simply "Popa") 
studied diazotrophically growing Anabaena oscillarioides filaments with a nanoSIMS 
(nanometer-scale secondary ion mass spectrometry) technique after they were provided 
with stable dinitrogen gas composed of ^^N. After four hours, they found dips of the fN 
distribution at the locations of heterocysts, and a noisy plateau between them (see red 
curve in Fig. [3b, below) — with no evidence of a concentration gradient consistent with 
diffusive transport away from heterocysts. 

It has not been immediately clear how to reconcile the Wolk and Popa fN 
distributions. It is tempting, but not necessary, to attribute quahtative differences 
between the Wolk and Popa distributions to distinct species-specific physiology. Instead, 
we hypothesize that fN transport is qualitatively similar in all heterocystous filamentous 
cyanobacterial species and investigate the quantitative fN patterns expected from 
modelling the two experimental approaches with a single model. 

Wolk et al [3] developed a quantitative model of fN transport and incorporation 
along the cyanobacterial filament, including diffusive transport and fN consumption 
in vegetative cells. However, the Wolk model did not distinguish periplasmic from 
cytoplasmic transport, used a homogeneous continuum model of the filaments, did not 
include stochasticity, and assumed that cytoplasmic fN was always rate-limiting for 
growth — i.e. growth vs local fN concentration was strictly linear. We found (see 
Appendix A) that with this last assumption the Wolk model was unable to reproduce 
the distribution of Popa after prolonged production of fN at heterocysts. More recently, 
a stochastic cellular computational model of cyanobacterial growth and fN transport 
was developed by Allard et al [22]. However, they assumed that only periplasmic 
fN transport was significant and they did not investigate the possibility of local fN 
starvation. We adapt their model, with both cytoplasmic and periplasmic transport 
but stripped of heterocyst development and supplemented by vegetative growth that is 
independent of non-zero fN concentrations, to the study of fN distributions. Using this 
model, we reconcile the two qualitatively distinct fN distributions of Wolk and of Popa. 
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Figure 1. Schematic of fixed-nitrogen dynamics as represented by Eqns. [1] and [2j 
Fixed- nitrogen amounts are indexed to each cell by i, with Nc{i) the cytoplasmic 
fixed- nitrogen and Np{i) the periplasmic fixed- nitrogen of each cell. Each shaded green 
region represents the cytoplasms of two adjacent cells. Dc describes the cytoplasm- 
cytoplasm transport, Dp the periplasm-periplasm transport, Di the import from the 
periplasm to the cytoplasm, De the export from the cytoplasm to the periplasm and 
Dl the loss from the periplasm to outside of the filament. 

We then use the quahtative Popa results to constrain possible periplasmic fN transport. 
We find periplasmic fN transport is not necessary to explain the current data but cannot 
be ruled out as a significant contribution. 

2. Model 

2.1. Fixed-nitrogen transport 

Our fixed-nitrogen (fN) transport and incorporation model is similar to that of Allard 
et al [2^, and tracks the total amount of cytoplasmic fN, Nc{i, t), in each cell i vs. time 
t: 

j^Ncit, t) = ^c{^) + DjNp{t, t) - DENcit, t) + G„ (1) 

where ^c{i) = ^Rcii — ^) + ^Lc{'i + ^) — ^Lc{i) —^Rcii) is the net diffusive cytoplasmic 
flux into cell i and is equal to the sum of the fluxes into the cell from the left and 
right minus the fluxes out of the cell to the left and the right. Similarly we track the 
corresponding periplasmic flxed-nitrogen Np{i,t): 

jNp{i, t) = $p(0 - {Di + Dl)Np{i, t) + DeNc{i, t), (2) 

where, similarly, $p(i) = ^Rp{i — 1) + $lp(^ + 1) — ^LP^i) — ^Rpii) is the net diffusive 
periplasmic flux into the periplasm of cell i. In addition to fluxes along the filament, Dj is 
the coefficient for periplasm-to-cytoplasm transport, De is the coefficient for cytoplasm- 
to-periplasm transport, and is the coefficient for losses from the periplasm to the 
external medium. These transport processes are shown schematically in Fig. [1] Note 
that any fN transferred between cytoplasm and periplasm, with Dj or De, shows up 
with opposite sign in Eqns. [T]and[2] — so that the total amount of free fN is unchanged 
by the exchange. In addition to these transport terms, there is a source/sink term Gi in 
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Eqn. [T] that describes fN production and consumption in the cytoplasms of heterocysts 
and vegetative cells, respectively. This G term is discussed more extensively below in 
Sec. [221 

Diffusive transport between cells is described by the flux terms $. Following AUard 
et al [22j and MuUineaux et al [TT] we assume that they are limited by the junctions 
between cells. Supporting this, there is no evidence in the images of fN by Popa et al [21] 
or of calcein by Mullineaux et al |TT] of subcellular gradients. Pick's law states that net 
diffusive fluxes between two cells will be proportional to the density difference between 
the cells. Each cell will have two outgoing fluxes, one to the left neighbour and 
one to the right neighbour each proportional to the local density N{i,t)/L{i,t) and 
to transport coefficients Dq ot Dp for cytoplasmic or periplasmic fluxes, respectively. 
Each cell will also have two incoming fluxes from its neighbouring cells. 

The value of the transport coefficients Dc or Dp will be proportional to the 
microscopic diffusivity of fN, but will also depend on the size and shape of cytoplasmic or 
periplasmic connections between cells — which are not well characterized. Nevertheless, 
we can take the cytoplasm-to-cytoplasm transport coefficients for calcein [11] and simply 
scale them to glutamine under the assumption that contributions due to connection 
geometry are unchanged. Calcein transport coefficients reported by Mullineaux et 
al [11] assumed a constant cell length. Multiplying their E in units of by a 
cell length of 3.38 fim, gives our transport coefficients in units of fim s~^. We 
then scale from calcein to glutamine with a factor of the cube-root of the ratio of 
molecular masses, assuming Stokes-Einstein diffusivity that inversely depends upon 
radius. We obtain D(7=1.54 /xm s~^ between two vegetative cells and Dc=0.19 fj,m 
s~^ between a vegetative cell and a heterocyst — we keep these values fixed in 
this paper. Note that with uniform transport coefficients, such as among clusters 
of vegetative cells, the flux terms in in Eqn. [1] or [2] look like a discrete Laplacian, 
VjiN/L) = N{i + 1, t)/L{i + l,t) + N{i - 1, t)/L{i - 1, t) - 2N{i, t)/L{i, t)) and lead to 



the more familiar diffusion equation when coarse-grained (see Eqn. [5] and Appendix A 
below) . 

We can estimate a lower bound of periplasmic fN transport with the GFP-FRAP 
data of Mariscal et al [15]. They observed a fluorescence recovery after photobleaching 
timescale of approximately 150s [I^. This corresponds to Dpgjp ^ L/150s ~ 0.02/im 
s~^, using a cell length L ^ 3/im. Scaling the transport of GPP (27 kDa) to glutamine 
(146 Da) by the cube-root of the molecular masses [23], i.e. by a factor of (27 kDa/146 
Da)^/'^ ^ 5.7, and taking this as a lower bound, then gives Dp > 0.1 /xm/s. However, 
the relatively large GPP may diffuse particularly slowly through periplasmic connections 
compared to glutamine. PRAP experiments with GPP in Anabaena PCC 7120 showed 
recovery in part of a single cell's periplasm after around 60s [16] and after around 5s 
^15j . This faster time of 5s is similar to the 4s recovery time observed in the periplasm 
of E. coli ^24j, which also measured an associated diffusivity of Dcoh ~ 2.6/im^/s. We 
use this diffusivity to determine an upper bound of Dp by scaling it to glutamine and 
dividing by cell length to obtain Dp = Dcou x 5.7/3.38/im ^ 4.4/im/s. The broad range 
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we consider for periplasmic transport is then 0.1 fim s^^ < Dp < 10 fim s~^. 

When we consider periplasmic transport, we must also consider exchange between 
the cytoplasm and the periplasm — as well as loss from the periplasm. We expect 
that uptake of fN from the periplasmic to the cytoplasmic compartment, controlled 
by Dj in Eqns. [H and [2], is an active process, controlled by e.g. ABC transporters 
[20] . Nevertheless, we expect that the rate of import is proportional to the periplasmic 
density, Np/L, and to the number of transporters, which will themselves be proportional 
to the cell length ~ L. Similarly, leakage of fN from the cytoplasmic to the periplasmic 
compartment, controlled by D^, will be proportional to the cytoplasmic density Nc/L 
but also to the amount of membrane or of (leaky) ABC transporters (~ L) [25]. The 
result is transport terms Dj and that are independent of length L and have units of 
s~^. Using the transport data of Pernil et al [20] for uptake of extracellular glutamine, 
we estimate (see Appendix B) Dj = 4.98 s~^ and De = 0.498 s~^. We also consider 
the possibility of stronger transport, with Dj = 49.8 s^^ with De = 4.98 s~\ because 
the outer membrane permeability barrier [17] may decrease extracellular uptake with 
respect to the periplasmic uptake Dj. There is evidence that fixed-nitrogen is lost from 
the filament to the external medium [261 EZ] and we include a loss term for the periplasm 
with coefficient D^ in the Eq. [2l We take the loss rate as very small compared to the 
reuptake rate, using 1% of the lower Dj value i.e. D^ = 0.0498 s^^ — qualitatively 
corresponding to the observation of a permeability barrier in the outer membrane |17j . 
Our results are qualitatively unchanged in the range of < Dl/Dj < 0.1. 

2.2. Cell Growth and Division 

We take PCC7120 cells to have a minimum size of Lmm = 2.25yum and a maximum 
size of Ljnax = '^Lmin [D IS]- When a cell reaches Lmax it is divided into two cells of 
equal length, each of which is randomly assigned a new growth rate and half of the 
fixed nitrogen in the parent cell. We found that average concentration profiles were 
largely insensitive to stochastic effects. Accordingly, most of our results came from a 
deterministic model where all cells were initialized with the same length Lq = 2.8/im 
and the same doubling time Td = 20h [191 EH]. For our stochastic results, presented 
in Figs. Et^b) and (c), we initialized lengths randomly from an analytical steady state 
distribution of cell lengths ranging between Lmin and L^ax [29]. We define a minimum 
doubling time T^m = Td — A and a maximum doubling time Tmax = Tp, + /S.. These 
then define a minimum growth rate Rmin = Lmin/Tmax and a maximum growth rate 
Rmax = Lmin/Tmin and we randomly and uniformly select optimal growth rates 
from the range Ri G [Rmin, Rmax]- The standard deviation of the maximal growth 
rate is then ap = a/2/3 ALmm/[(^z) + ^)(^d — A)]- We expect that the coefficient of 
variation of the growth rate ap/Ravg is of similar magnitude to that seen in a study of 
mutant Anabaena [22], which was crp/Ravg ~ 0.165. 

In Eqn. [H for vegetative cells Gi = G^eg is a growth term determining the rate of 
freely diffusing cytoplasmic fN removed to support growth of the cell. Gyeg depends on 
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the actual growth rate R of the cell (in fim/s), which in turn depends upon the locally 
available cytoplasmic fN: 

G,eg = -gR{R°''\Ncii,t)), (3) 

where g is the amount of fixed nitrogen a cell needs to grow per /im - which we now 
estimate. Dunn and Wolk [30] measure a dry mass for one cell of A. cylindrica of 
1.65xl0~^^g. Cobb et al [31] measured the nitrogen content of A. cylindrica to be 
5-10% of dry mass (see also [32]). We take their typical value of 10% [31] and scale 
it by volume from A. cylindrica to PCC7120. A. cylindrica appears to be 1.5 x as 
wide and 1.5x as long as PCC7120 [1], giving ^ 2.07x10^° N atoms per average cell 
of PCC7120. According to Powell's steady state length distribution of growing bacteria 
[29] the average cell size is 1.44x as long as the smallest cell, so that ^ 1.4x10^^^ 
N atoms are needed for a newly born cell to double in length. This implies that 
c/ = 1.4 X lO^y Lmin ^ 6.2 X lOVm^^ 

We assume that cells increase their length at their optimal growth rate as long 
as there is available cytoplasmic fixed- nitrogen (see discussion in Sec. Hj). Otherwise, 
they can only grow using the fixed-nitrogen flux into the cell from neighbouring cells 
and the periplasm: 

^jRr\ ifiVc(M)>o 

\ min{^,Jg, Rf) if Nc{i, t) = 

Note that cells with Nq = may still grow, but will be limited by the incoming fluxes 
of fixed-nitrogen from adjoining cells and the periplasm, G^eg = ^in = ^Rc{i ~ 1) + 
^Lci''' + 1) + DiNp{i). When the incoming flux exceeds the requirement of maximal 
growth, then the cytoplasmic nitrogen concentration will rapidly become non-zero since 
we always have R < R"^^ . 

In Eqn. [H for heterocyst cells Gi = Ghet- The heterocyst fN production rate Ghet 
is chosen to supply the growth of approximately 20 vegetative cells and is 3.15x10^ 
unless otherwise stated. 



2.3. Some numerical details 

Filaments were initiated with two heterocysts, separated on either side by vegetative 
cells, on a periodic loop. Periodic boundary conditions were used to minimize end 
effects. Vegetative cells were then allowed to grow and divide, subject to available fixed- 
nitrogen. Because of strong spatial gradients of free (unincorporated) fixed-nitrogen 
between heterocysts, we group and report our data with respect to the number of 
vegetative cells between two heterocysts — typically we took a separation of 20 cells. 
This is approximately twice the typical heterocyst separation [Il[2], which allows us to 
investigate fixed-nitrogen depletion at the midpoint — corresponding to the location of 
an intercalating heterocyst. 

We distinguish free from incorporated fixed-nitrogen (fN). Free fN is freely diffusing 
in the cytoplasm, and is simply given by Nc{i,t). We report this as a linear density 
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Pp = Nc{i)/Li. Incorporated fN is the fixed-nitrogen that has been incorporated by 
cellular growth through the growth term G^eg- We report the incorporated concentration 
Pi = f Gvegdt/Li, where the growth is integrated over a fixed duration. During cell 
division, half of any previously incorporated fN is assigned to either daughter cell. Note 
that since we are interested in isotopically labeled fixed-nitrogen [3], El], we set all 
Nc{i) = at the start of our data gathering (t = 0), corresponding to the introduction 
of isotopically labeled dinitrogen to the filament that is then fixed by the heterocysts 
and subsequently supplied to the filament as free fixed-nitrogen via Ghet- The total 
fN, for e.g. Fig. HJ in a cell is p^ = Pf + Np/L + pj, where we include periplasmic 
contributions as well. 

3. Results and Discussion 

3.1. Short-time fixed-nitrogen distributions 

We first examined systems with only cytoplasmic transport (i.e. Dj = Dp = Dp = 0) 
at short times comparable to those of Wolk et al |i3j . In Fig. [2] we show fixed-nitrogen 
(fN) distributions along the filament, with respect to distance in cells from the nearest 
heterocyst, i. The cytoplasmic transport parameter was held constant at Dc = 1.54 
/im/s, estimated from calcein transport values in Section [2711 We varied the time-interval 
between the start of labeled nitrogen fixation {t = 0) and the times {t = 10s, 30s, 60s, 
120s, and 180s — see legend in panel (6)) at which the fN concentration profile along 
the filament was considered. As shown in (a), the incorporated fN, p/, increases steadily 
over time for all vegetative cells reached by free fN, while incorporated fN remains at 
zero in the non-growing heterocyst ("H"). As shown in (b), the free fN, pp, spikes at the 
heterocyst, where it is produced, and progressively decreases with distance away from 
the heterocyst. At later times, the free fN approaches a steady-state — the distribution 
does not change significantly with time. We show in (c) the total fN, pr = Pi + Pf- At 
very short times the free fN dominates the total, resulting in the spike in the distribution 
at the heterocyst. At longer times the incorporated fN dominates the total — which 
then loses its spike at the heterocyst and instead begins to form a dip at the heterocyst 
(see below). 

The radiographic technique of Wolk et al [3] does not distinguish between free and 
incorporated fN, so it is appropriate to compare it with the total fN pp from Fig. [2](c). 
In (d) we show the average of the two distributions determined by Wolk et al p], with 
filled red circles (left scale). We superimpose a qualitatively similar pp curve at t = 30s 
with the blue curve (right scale). We have shifted the pp = axis to agree with the 
experimental background far from the heterocyst, and have adjusted the scale so that the 
curves agree at the heterocyst. The agreement is not bad, and captures the qualitative 
feature of the Wolk data — namely the distinctive spike at the heterocyst. This spike is 
a feature of short exposure to labeled fN. Indeed, while Wolk et al took data at t = 120s 
we find our best qualitative agreement at t = 30s. In our model heterocysts immediately 
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Figure 2. Model fixed-nitrogen (fN) concentration vs. cell index i (the number of 
cells from a heterocyst at ? = 0), for short intervals of labeled nitrogen fixation ranging 
from 10s to 180s as indicated by the legend in (b). The incorporated fN pj is shown in 
(a), the free fN pp in (b), and the total fN pT in (c). In (d) the red circles indicate the 
experimental data digitized from Wolk et al [3] , and the blue squares are the t — 30s px 
data from (c) (using the right axis, scaled and shifted to best agree with experimental 
data). Our model data is from our deterministic model with no periplasmic transport, 
i.e. Dc = 1.54 p.m s'^ and Dj = Dp = Dp = 0. 



provide fN able to be transported. This is not the case in real heterocysts: glutamine 
levels in Anabaena cylindrica linearly increase during the first minute of exposure to N2, 
and other fixed nitrogen products appear after two minutes [6]. This one to two minute 
lag between exposure to labeled fN and its transport and incorporation is consistent 
with the 90s difference between the illustrative curves in Fig. EJ^d). 

3.2. Steady- state fixed-nitrogen distributions 

As seen in the last section, for t > 120s the free fN distribution has reached an 
approximate steady-state. We can obtain an analytic expression for this steady-state 
distribution of pi? by approximating the filament as a continuous diffusive medium. The 
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diffusion and consumption of free fN witliin a continuous filament is 

dpp/dt = DV^pf - c (5) 

witfi pf as the free fN concentration; D is the diffusivity and c is the fN consumed 
by growth per unit length. For the steady state {dpF/dt = 0) the general solution is 
quadratic (parabolic) vs. distance x from a nearby heterocyst. In steady-state, the 
flux on either side of the heterocyst at x = equals half of the fN production, i.e. 
Ghetf^ = —DdpF/dx. This flux, Ghet/'^i must also equal the total fN consumed by the 
growing vegetative cells on either side of it, i.e. cxq where xq is where growth halts 
and Pf{,xq) = 0. Applying pFiGhet/C^c)) = 0, and -DdpF/dx{0) = Ghetf^ gives us the 
steady-state profile 

Pf{x) = X X -\ — , (6) 

'^^^ ' 2D 2D 8cD' ^ ^ 

for X < Xo and Pf{x) = for x > xq. This is in excellent agreement with the steady- 
state Pf from our numerical model with non-stochastic cytoplasmic-only transport, as 
shown in Figure [3]^a). 

In our non-stochastic model all cells with free fN levels above zero grow at the same 
rate and so incorporate fN at the same rate. This produces a constant plateau in pj for 
vegetative cells between the heterocyst at x = and the cell where the free fN vanishes, 
i.e. for < X < Xq. This is seen in our discrete cellular model, as shown with px in 
Fig. [3](a). We note that the cells at i = 9 and i = 11 grow at a reduced rate despite 
having pf = 0, since some free fN is provided from their neighbours — i.e. > in 
Eq. m The plateau height depends upon how much time has elapsed since labeled fN was 
introduced. In Fig. [3](a) we show results for t = 4h, to be able to directly compare with 
the data of Popa [21]. The predicted plateau height {<R"P^>gt / <L>) from our previous 
continuum model is shown with a black dashed line, and the agreement is excellent. 
We also see that the free fixed-nitrogen pf (scale on right side) is approximately two 
orders of magnitude smaller than the total fixed-nitrogen px (scale on left side), which is 
dominated by the incorporated fixed- nitrogen pj. We see excellent qualitative agreement 
with the plateau-like distribution seen by Popa. The expected gradient of free fixed- 
nitrogen is too small to be easily resolved by their technique, compared to the much 
larger plateau of incorporated fixed-nitrogen. 



3.3. Stochasticity of Long Time Distributions 

Fig. El^b) shows the total fN distribution using our stochastic model that has been 
averaged over 1000 vegetative segments between two heterocysts. Qualitatively it looks 
like our deterministic model results. In (c) we show both the experimental results 
of Popa et al [21] (red curve, left axis) and one single sample distribution from our 
stochastic simulation (blue curve, right axis, with an arbitrary factor to separate the 
curves). Qualitatively the distributions are very similar, but strikingly rougher than the 
smooth curves in (a) and (b). The pf in (c) is smooth and relatively small in magnitude, 
indicating that the variation is almost entirely due to incorporated fN. 
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Figure 3. (a) non-stochastic concentrations after t — 4h of labeled fixed-nitrogen (fN) 
exposure vs. the cell index i. Heterocysts are indicated by ("H"). Shown are the total 
fN {pT, blue solid line with circles using the left axis) and the free fN (pF, red solid line 
with squares using the right axis). Note the dramatic difference in scale. The steady- 
state analytic result for the free fN from a continuum approximation, Eqn.[6l is shown 
as a green line; while the corresponding plateau in px due to incorporated fN is shown 
by the black dashed line, (b) averaged stochastic pr using doubling-time variation A 
= 4.5h, with error bars too small to be seen outside data points. The average is over 
1000 distributions, (c) the experimental fN concentrations of Popa et al [2T] in red 
circles (scale on left axis) and from a single stochastic result from our model using A 
= 4.5h and Ghet = 4.25x 10^ s^^ with px in blue squares and 25pF in green diamonds 
(scale on right axis, offset for clarity). Also shown as black lines are the vegetative cell 
averages, (d) the coefficient of variation in vegetative cell fN levels, ar/^Pr, vs. the 
coefficient of variation in the growth rate of the vegetative cells, au/Ravg- The error 
bars indicate statistical errors of the average variation measured over 100 samples. The 
grey shaded region indicates the size of the standard deviation, which corresponds to 
the expected error for a single sample. The red dashed line indicates the variation in 
the experimental data of Popa et al [21], with a single sample. For (a)-(d), we have 
used only cytoplasmic transport with Dq = 1.54 /im s""'^ and Dj — De ~ Dp = 0. 
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In our stochastic simulation, the variation of pt shown in (c) is due both to 
the varying growth rates of cells along the filament and to their varying lengths — 
which leads to corresponding differences in the density of incorporated fN after a fixed 
time. This raises the possibility that the variability seen by Popa et al [21] is not due 
to instrumental or sampling noise, but simply reflects growth rate variability in the 
filament. In (d) we plot the coefficient of variation of the total-fixed nitrogen, ax/^Pr, 
vs the coefficient of growth rate variation, crR/Ravg- Here aji is the standard deviation 
of the growth rate distribution, and Ravg is the average growth rate. Similarly ctt is 
the standard deviation of pt in vegetative cells, while Ap^p is the difference between 
the average pr in vegetative cells and in heterocysts. While an/ Ravg was calculated 
exactly, we measured otI Ap^ over 100 segments of cells, and used this to determine 
the average coefficient of variance, the statistical errors, and the standard deviation 
around the average — as shown in (d). With the dashed red line, we also show ar/ Apr 
from the Popa experimental data in (c). The grey region shows one standard deviation, 
which approximates the single interval error for the Popa data. We see that a broad 
range of growth rate stochasticities are consistent with the Popa's observed variation, 
including the coefficient of variation ctr/ Ravg = 0.165 reported for a mutant of PCC7120 
[22] . We conclude that the variability seen by Popa et al is probably determined by the 
stochasticity of the growth rate. 

After extended isotopic dinitrogen exposure, much longer than typical cell-doubling 
times, we expect a nearly constant distribution of pr — corresponding to the nearly 
uniform concentration of fN in a cyanohacterial filament. Notably, the variation within 
the plateau, evident in the data of Popa et al [21] and modelled in Fig. [3](d), will also 
vanish at late times — since it arose due to differential growth rates leading to variations 
in labeled vs. unlabelled p/. Only old heterocysts, which retain some unlabelled fN, will 
present local dips in the concentration profile. On top of this nearly constant p/, will be 
much smaller steady-state free fN gradients coming from heterocysts — as described by 
Eq. [6]and illustrated in Fig. [3]^a-c). Imaging fN distributions in filaments that have been 
grown for weeks in labeled dinitrogen may allow small gradients of pi? to be discerned 
with NanoSIMS techniques [2^. 

3.4- Long Time Nitrogen Distributions With Periplasmic Transport 

We have established that the qualitative characteristics of both the spiked of Wolk et 
al [3] and the fiat distribution of Popa et al [21] can be reproduced using a single 
quantitative model for the cyanohacterial filament using only cytoplasmic transport. 
In this section, we explore the question of how much additional periplasmic transport 
is possible while retaining the qualitative phenomenology. The strongest effect is seen 
after a Popa-like plateau is formed in p^. 

In Figs, m^a) and (b) import and export between the cytoplasm and periplasm has 
Di = 4.98 s~^ and De = 0.498 s~^ (see Appendix B). The only parameter that varies 
is the value of the periplasmic diffusivity Dp, as indicated by the legend in (b). We 
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Figure 4. Fixed-nitrogen distributions for systems with additional periplasmic 
transport, as discussed in Section 13.41 In all four plots the periplasmic transport 
rate Dp (in units of fim s^^) is varied according to the legend in (b). (a) the total 
fixed-nitrogen distribution, pT, for transport between the cytoplasm and periplasm 
with Dj = 4.98 s~^ and De = 0.498 s~^. The data of (b) corresponds to the curves 
of the same colour in (a) and shows the fraction of the total fixed-nitrogen fiux along 
the filament that is periplasmic vs. the cell index i. (c) and (d) are similar to (a) and 
(b), respectively, except that import and export between the cytoplasm and periplasm 
are tenfold stronger, with Dj = 49.8 s~^ and De = 4.98 s~^. As before, £'c=l-54 fim 
s-\ We use Dl = 0.0498s-i and Ghet = 3.67xl0'3 s-\ 

vary Dp up to the magnitude expected for unimpeded diffusive transport of glutamine 
along the filament, as discussed in Sec. 12.11 The characteristic dip in the centre of 
the total fixed-nitrogen distribution, px, in (a) fills in at larger values of Dp — but 
maintains its qualitative features throughout the range. The mid-segment dip fills in 
because periplasmic transport allows free fixed-nitrogen to move past the first starving 
vegetative cell with Nc = 0. To evaluate how significant the periplasmic transport 
is, we show in (b) the fraction of the flux along the periplasm to the total fiux along 
the cytoplasm and periplasm. The fiux is mostly cytoplasmic next to the heterocysts, 
and mostly periplasmic next to starving cells. In between, up to 50% of the fiux can 
be periplasmic. For the larger Dj and De values of Figs. ID^c) and (d) the results are 
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very similar. We conclude that significant, though probably not dominant, periplasmic 
transport is consistent with the qualitative fixed-nitrogen distribution reported by Popa 
et al [21]. 

4. Summary 

We have presented a quantitative model for fixed-nitrogen (fN) transport and dynamics 
in a cyanobacterial filament, including growth of vegetative cells and production of fN 
by heterocysts. Our model reproduces the qualitative fN distribution near heterocysts 
(Fig. |2]) seen after short exposure to isotopically labeled dinitrogen by Wolk et al. 
[3] ('Wolk"); but also reproduces the qualitative distribution between heterocysts 
(Fig. [3]) seen with nanoSIMS techniques after much longer exposure by Popa et al. 
pTj ( "Popa" ) . The apparent lack of a fN gradient in the Popa distribution is explained 
by the much larger amount of incorporated fixed-nitrogen, which gives the characteristic 
noisy plateau between heterocysts. Although much smaller in magnitude, the modelled 
free fN distribution does show a smooth gradient away from heterocysts (Fig. [3]^a)) 
refiecting diffusive transport. 

Qualitatively, our results indicate the importance of the large, noisy, plateau- 
like concentration profile of incorporated fixed-nitrogen. The level of the plateau will 
increase linearly in time after the free fN gradient away from heterocysts reaches its 
steady-state, which takes minutes. The free fN, pp, is soon much smaller than the 
incorporated fN, p/. The spike in the Wolk distribution, seen in the first few minutes 
of exposure to isotopically labeled dinitrogen, was when the free fN still dominated 
the total pt- In contrast, the plateau seen in the Popa distribution was after hours of 
exposure, when the total fN was dominated by the incorporated component and the 
gradient va. pp was relatively small even with respect to variations of pi within the 
plateau (see Fig. Et^c)). At these intermediate times, long after the steady-state pp but 
much shorter than typical division times, the heterocysts present dips of labeled fixed- 
nitrogen — as seen by Popa [21] . Stochasticity, through initial cell lengths and cellular 
growth rates, adds noise to the incorporated fN distribution, p/. As we show in Fig. [3] 
(d), the modelled cell-to-cell variation is due to differences in growth rate which manifest 
themselves as different amounts of labeled fN incorporated by growth. We believe this 
explains the observed variation in the experimental nitrogen distribution of Popa et al 

The possibility of periplasmic fixed-nitrogen transport raised by a contiguous 
periplasm [TJ], transport of GFP [T5], and an outer membrane permeability barrier 
|17j . led us to explore the possible impact of periplasmic transport in addition to 
cytoplasmic cell-to-cell connections [H]. As illustrated in Fig. HI we found that a 
significant proportion (up to 50%) of fN transport along the filament can be observed 
without qualitatively changing our fN distributions. We cannot rule out a role for the 
periplasm in fN transport, and larger proportions may be achievable with more extensive 
parameter searches. Indeed, the fact that knocking out amino-acid permeases leads to 



Reconciling cyanohacterial metabolite and transport experiments with modelling 



15 



impaired diazotrophic growth IT^ EU] does indicate that significant leakage from 
the cytoplasm does occur (i.e. De > 0) — which is a necessary part of periplasmic 
transport. Nevertheless, our results are also consistent with no significant periplasmic 
transport. 

We have used a simple model, Eq. HI to limit growth in response to local fixed- 
nitrogen starvation: cells grow at a fixed rate R"^* if there is sufficient free fixed-nitrogen 
available locally [Nq > 0), and at the maximal rate allowed by local input fluxes of free 
fixed-nitrogen if not. This model reflects the absolute limitation on growth placed 
by available fixed-nitrogen, but also recovers the characteristic plateau of pt seen in 
the Popa distribution. As we illustrate in Appendix A "graded" growth rates, that 
smoothly depend on Nq and vanish at Nq = 0, respond to the graded distribution of 
free fixed-nitrogen away from heterocysts, Eq. IA.21 and produce graded growth patterns 
and graded distributions of pi and pt (see Fig. lAll) — unlike those observed by Popa [2T] . 
Graded growth models, such as the Monod model [33] , are typically designed to describe 
response to external metabolite concentrations rather than cytoplasmic concentrations. 
Like our growth model. Salmonella typhimurium [31] and Escherichia coli [351 ES] both 
show non-graded growth in response to cytoplasmic fixed-nitrogen starvation. However, 
there is little direct evidence of how growth of vegetative cells in cyanobacterial filaments 
depends on the cytoplasmic freely available fixed-nitrogen, pp. 

The experimental work was done on different species of filamentous cyanobacteria, 
namely Anabaena cylindrica by Wolk et al [3| and Anabaena oscillarioides by Popa et al 
[2T] . In contrast, our models were parameterized by consideration of transport studies 
done in Anabaena sp. strain PCC7120 (see Sec. 12. ip . While we expect parameter 
values to differ between species, so that precise numerical agreement with the Wolk 
and Popa results is not to be expected, our qualitative agreement indicates that similar 
fixed-nitrogen transport and dynamics could apply in these different model organisms. 

We expect the qualitative features of the Wolk and Popa fixed-nitrogen distributions 
to apply to the early and late-time fixed-nitrogen distributions, respectively, of all 
filamentous cyanobacteria with heterocysts that are growing in media lacking fixed- 
nitrogen and that have significant cytoplasmic connections. Nevertheless, quantitative 
details will depend, via the model parameterization, on the cyanobacterial strain and 
the experimental conditions. For example, sufficiently large periplasmic transport (via 
Dp) would broaden Popa's dip in incorporated fixed-nitrogen between heterocysts due 
to transport and growth past the first starving cell. Conversely, a larger cytoplasmic 
diffusivity (via Dc) would alter the magnitude of Wolk's free fixed-nitrogen (via Eqn. [6]) 
but not the incorporated fixed-nitrogen. The average doubling time determines the 
characteristic time at which we expect the later time distribution to emerge. Cell-to- 
cell variations in the doubling time, in turn, determine the amount of variation expected 
within the Popa-like plateau of incorporated fixed-nitrogen near heterocysts for the first 
few generations of labelled growth. 

Our results reconcile the metabolite patterns shown by the autoradiographic 
technique of Wolk et al [3] with the nanoSIMS technique of Popa et al [21j. It is 
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not obvious which would currently have the best resolution, and it would be useful to 
have either or both applied to Anabaena PCC7120 — the most popular current model 
for filamentous cyanobacteria and for which recent transport studies exist [HI [15]. We 
can make four qualitative but testable predictions from our work. First, as the duration 
of labelled fN is increased the total fN distribution away from heterocysts will evolve 
from a Wolk-like peak, dominated by free fN, to a Popa-like plateau, dominated by 
incorporated fN. As indicated by Fig. |2l the crossover takes place after only a few 
minutes of labeled dinitrogen fixation. Second, on a similar timescale the free fN will 
reach a steady-state smooth parabolic shape even in a single sample. This free fN may 
be directly measurable, for example, by fluorescence resonance energy transfer (FRET)- 
based metabolite nanosensors [37]. The parabolic gradient of pp could also be uncovered 
by averaging the results of many experiments measuring pT- Third, the variations seen 
by Popa et al [21] will decrease in relative magnitude after at least one doubling-time 
of exposure to labeled fN. Fourth is more of an assumption, but we expect that growth 
of vegetative cells will slow abruptly when the freely available fN vanishes — according 
to Eqn. m rather than in a graded fashion according to Eqn. IA.3I Finally, with more 
precise parameterization of our transport parameters, sources, and sinks we can make 
quantitatively testable predictions. Prediction of the variance of cell-to-cell fN, aj, 
requires only the observed growth rate variance and the duration of exposure to 
labeled dinitrogen. 
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Appendix A. Graded growth models 



In this Appendix we address the effects of a graded growth model, where the growth rate 
of a cell depends continuously on the cytoplasmic free fixed-nitrogen concentration pp = 
Nc/L, and vanishes when Nc = 0. Examples of such graded growth models include the 
Monod model (normally applied to concentrations of extracellular metabolites) [33] and 
the growth model used by Wolk et al [3j for modelling the response of cyanohacterial 
filaments to short exposures to labeled dinitrogen. Here we consider the model of Wolk 
et al [3], though we expect similar results with any graded growth model. 

The graded-growth model of Wolk et al [3j has growth linearly proportional to free 
fixed-nitrogen concentration, pp. The dynamics of pp are then 

dpp/dt = DV'^pf - kpF, (A.l) 

where the first term on the right represents diffusion along the filament and the second 
term is the local sink due to vegetative growth. (Note the difference with our coarse- 
grained growth model in Eq. [H) We ignore the discrete cellular structure of the filament 
in order to obtain an analytic solution. A heterocyst at x = will impose a flux 
J = —Ddpp/dx = Ghetf^ there, leading to a steady-state solution 

PF = :^^e-v^-. (A.2) 



2VkD 

The incorporated fixed-nitrogen is proportional to this, and so pr will exponentially 
decay away from the heterocyst with a length-scale ^JUJk. This is not observed in the 
distribution of Popa et al [21] — since the plateau shown in Fig. [3] is inconsistent with a 
small \/ Djk while the sharp dip seen between existing heterocysts is inconsistent with 



a large ^JD jk. 
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Figure Al. The numerical total fixed-nitrogen distributions vs. cell index i using 
graded growth proportional to local pp as Eq. IA.31 The growth-rate proportionality 
constant k (in units of \xva. s""*^) is varied, as indicated by the legend. This is a 
deterministic cellular model with cytoplasmic transport: Dq — 1-54 /xm s~^ and Di 

^ De = Dp = 
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With our discrete deterministic model, we have modelled the Wolk-style growth 

RgradediN,) = kNjU (A. 3) 

The results are shown in Fig. lAll for various values of k. The discrete nature of the 
model modifies the observed pattern near the bounding heterocysts compared to Eq. IA.2I 
— but either a plateau but no mid-segment dip is recovered (for small k) or a wide mid- 
segment trough with no plateau is recovered (for large k). We are unable to recover the 
experimental distribution seen by Popa et al ^1] with such a graded growth model. 

Appendix B. Parameter estimation of Dj and De 

To estimate the value of Dj (from Eqns. [T] and [2]) we begin with the measurement of 
Pernil et al |20j that 1.65 nmol (mg protein)"^ is imported into a PCC7120 filament 
in 10 minutes from a medium containing 1 /xM glutamine. From our discussion in 
Section 12.21 the dry mass of a PCC7120 cell is 4.89xl0~^^g and the mass of protein, 
approximated at 55% of the dry mass of the cell [38], is 2.69xl0-^2g Tj^jg implies that 
there is 4.44x10"^ nmol of glutamine, or Nq = 2.67x10^ glutamine molecules, imported 
in 10 minutes into each cell. The term in Eqn. [1] describing import from the periplasm 
to the cytoplasm is dNc/dt=DjNp. If we assume that the periplasmic concentration is 
equal to the glutamine concentration in the external medium, then Np can be replaced 
by the product of the glutamine concentration in the periplasm, pp, the cross sectional 
area of the periplasm A (for a O.l/xm thick periplasm surrounding a cytoplasm 1/im 
in radius A = 0.2l7ryum^), and the cell length L giving Nq = DjppLAr, with r = 10 
minutes. Solving this for Dj using L = Lmm = 2.25/xm yields Dj = 4.98 s^^. ABC 
transporters, known to transport amino acids [19], are asymmetric, and we assume that 
export is lOx weaker than import, and that De = 0.498 s~^. It is also known the the 
outer cell membrane provides a permeability barrier [T7], implying that the periplasmic 
glutamine concentration may not be as high as the external glutamine concentration, 
and so we also consider a. De and Dj pair that is ten times larger with Di = 49.8 s~^ 
and De = 4.98 s-\ 
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